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Abstract 

In recent years a very exciting and intense activity has been devoted to the un- 
derstanding and construction of materials that enjoy exotic optical properties, such 
as a negative refractive index. Motivated by these experimental and theoretical 
developments, we use the string- inspired idea of holography to study the electro- 
magnetic response of a certain class of media: strongly coupled relativistic systems 
that admit a dual gravitational description. Our results indicate that this type 
of media generally have a negative refractive index. Moreover we observe that a 
negative refractive index could be a common feature of relativistic hydrodynamic 
systems at low frequencies. 
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Introduction 

It is a well-known fact that there exists a close analogy between the free propagation of 
light in curved spacetime and the propagation of light inside a material in flat space. The 
geometry bends the light trajectories along null geodesies of the spacetime metric gfj,^, 
while the electromagnetic properties of the medium deviate the trajectory of the light 
according to its electric and magnetic tensors: eij,fiij. Indeed it is formally possible to 
rewrite the microscopic Maxwell equations in curved spacetime as macroscopic Maxwell 
equations in a continuos medium in flat spacetime (see for example [1]). 

This old story takes a new turn with the recent discovery that one can engineer "real" 
materials with highly unusual and astonishing electromagnetic properties [2ll3]. They are 
called metamaterials, since they are artiflcial materials engineered to provide properties 
which are not readily available in nature. They are typically constructed using periodic 
arrays of very small and closely spaced elements, playing the role of atoms in natural 
materials. One can engineer materials that force the light to follow almost any desired 
trajectory, allowing the construction of cloaking devices, perfect lenses, photonic black 
holes and so on. We cannot do justice to the literature here, we refer the reader to the 
reviews [UllHT]. 

One of the most interesting properties of some metamaterials is the negative refractive 
indeJ^ The theoretical possibility of this occurrence was flrst envisaged in 1968 by Vese- 
lago [8], although there are earlier suggestions due to Lamb, Schuster, and Mandelstam, 
see [9] and [10] for a historical account. In such a material the phase of a wave packet 
propagates in the opposite direction to the energy flux, which in turn leads to a number 
of surprising phenomena: many familiar laws of optics are modifled, e.g. the Snell's law. 



^Materials with this property are also known as Left-Handed Materials. 
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the Doppler effect, the Cherenkov effect etc. Around 2000 the first samples of negative 
refractive materials (NRM) were built [2||3|[TT| and since then they have attracted intense 
interest, due to their multiple technological applications. 

It would be important to find new and general classes of physical systems that exhibit 
negative refraction, and for which it is possible to exactly compute some relevant physical 
quantities, at least at the linear response level. In this paper we observe that the negative 
refraction seems to be a generic prediction of relativistic hydrodynamics. More precisely, 
we show that if the correlator of the transverse current coupling to the EM field is domi- 
nated by a diffusive pole then there is a low-frequency region where the negative refraction 
occurs. The presence of the diffusive pole is a consequence of the hydrodynamic equations 
in the case of relativistic field theories at finite temperature and chemical potential. 

We study in detail the refractive index for the specific example of a strongly coupled 
medium (eventually coupled to a dynamical electromagnetic field) that admits a gravity 
dual weakly coupled description. The computation of the transport coefficients and of the 
thermodynamical quantities is performed by using the string theory-inspired approach of 
holography. The holographic analysis goes well beyond the hydrodynamic regime, because 
the response functions are calculated at every frequency. Our results show that certain 
strongly coupled systems, coupled to a dynamical electromagnetic field, have negative re- 
fractive index in some specified range of frequencies. This confirms the hydrodynamic 
argument and introduces a new connection between gravity and optics, based on the idea 
of holography, realized in string theory by the AdS / CFT correspondence [I2] . 

The AdS/CFT correspondence is a map between a quantum field theory in fiat space, 
and a classical field theory in a completely different spacetime - a curved one, and with a 
larger number of dimensions. Often this is an Anti-de Sitter space or some deformation 
of it. The correspondence is holographic, since one can interpret the fiat spacetime of 
the QFT as being the boundary of the curved spacetime. In this setup the geometric 
properties of the higher dimensional curved spacetime are translated in properties of 
matter in fiat spacetime (see [13] for a review). 

In the last decade, this holographic duality has been used with the aim to model 
several systems realized in the laboratories. The first example was the quark-gluon plasma 
produced in heavy- ion collisions [H]. This was followed by holographic realizations of 
many phenomena of condensed matter systems: superconductivity and superfluidity [15] , 
Fermi gas at unitarity [I6], the quantum Hall effect [17], non- Fermi liquids [18], quantum 
phase transitions [19j. See ^20j for reviews. 

This paper is a first attempt towards an optics/geometry holographic duahty, where 
string theory can provide a useful description for some unusual optical properties of elec- 
tromagnetic media. 

We believe that our results could be interesting both for the optics community, be- 
cause they offer new theoretical examples and a "laboratory" for negative refractive index 
materials, in which the electromagnetic properties of the medium can be computed ex- 
phcitly; and for the string theorists, because they provide a new exciting field of research 
in the AdS/CFT correspondence. 

We should mention at the outset a generic caveat of this string theory setup. The 
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systems we are studying are not automatically coupled to a dynamical photon, so there 
is strictly speaking no propagation of light; nevertheless we speak of the refraction index. 
We can do it because, at the first order in the electromagnetic coupling, the optical 
properties of the medium are completely determined by the linear response to an external 
electromagnetic field, which can just be taken as a background, non-dynamical field. At 
the phenomenological level we can think of our system as a strongly interacting medium 
weakly coupled to a dynamical photon. The first part of the system is studied using the 
holography, while the coupling to the photon is taken into account perturbatively at the 
linearized level. 

We should also notice that the systems we study are relativistic, homogeneous, isotropic, 
strongly-coupled media. We cannot pretend that our results can be directly applied to 
the real metamaterials that in particular are non-relativistic. It would be interesting to 
understand if our results can give experimentally testable predictions for systems that 
occur in nature (perhaps in neutron stars or other exotic situations). Moreover it could 
be interesting to study linear the EM response in non-homogeneous, non-isotropic or 
non-relativistic holographic setups. 

The paper is organized as follows: in section [1] we recall the formalism of electro- 
magnetic fields in a medium and the linear response theory that allows to derive the EM 
response tensors. In section |2] we discuss the conditions under which negative refraction of 
light occurs, and we give our main result: negative refraction is generic in a certain class 
of relativistic systems. In section |3] we describe the holographic setup and the techniques 
of the computation. In section H] we give more details on the results, then we conclude 
with some additional comments and future directions of investigation. In the Appendix 
we discuss some details regarding both the optics and the holographic setup. 

1 Macroscopic EM and linear response theory 

The analysis of the electrodynamics of continuos media is traditionally performed by 
introducing, in addition to the electric and magnetic field, E and B, the macroscopic 
fields D and H. They are not independent fields but contain the information about the 
response of the medium to the application of an external field. The general relation, at 
the linearized level, i.e. for weak external fields, has the form 



In this approach the properties of the medium are encoded in the (Fourier-transformed) 
response functions eij{u), Hijiu). They are often assumed to be frequency- dependent only. 
In an isotropic medium, moreover, they must reduce to scalar functions: eijiu) = e(w) 6ij, 

If the medium has spatial dispersion the response functions also depend on the wave 
vector k. In this case it is possible [21] to reduce the number of macroscopic fields to D, 




(1) 
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E, B, and use the single function eij{oj, k) to describe the hnear response of the medium: 

Di = eij{u,k)Ej (2) 

For isotropic media the dielectric tensor e^j can be decomposed in its transverse and 
longitudinal part as 

eij = eriuj, k) (^ij - + ei(u;, k)^^ (3) 

The dispersion relations are obtained by projecting the Maxwell equations on the trans- 
verse and longitudinal components. We have 

eT{uj,k) = —, eLiuj,k) = (4) 

In experimental optics the most common description is not in term of er and but it is 
in terms of e and fi. The two approaches can be related in the following way. At /c = 
one has = = €{u). If we expand the transverse dielectric permittivity as a series in 
k, we can write [9] 

eriu, k) = e{uj) + —Al- ^1 + 0{k') . (5) 



From dl]) and (E]), keeping the terms up to k"^, we find that the transverse dispersion 
relation can be written as 

k"^ 

— =n'^{uj) = e{uj)fi{uj) (6) 

where n{u) represents the refractive index, and this is the same relation usually obtained 
in the e — /i approach. The two approaches are equivalent when the spatial dispersion is 
small, and one can neglect the higher-order terms in the expansion ([5]). 

It is important to note that the n{uj) defined via the eq. ([5]) is not the same that 
appears in ([T]). It is instead an effective containing both electric and magnetic effects. 
The magnetic permeability, as commonly defined, can be obtained from the response 
functions as |21] 

1 1 2 r eT{uj,k) - eL{oJ,k) 

1 = ^ limfc^o n • (^) 

fi[ujj k^ 

However it is the effective fi that correctly reproduces the dispersion relation and hence 
the electromagnetic propagation and dissipation!^. 

The electric permittivity and the magnetic permeability can be obtained from linear 
response theory [22]. In the linear response theory the electromagnetic current Jj is 
proportional to the vector potential Aj, Ji = GijAj, where Gij is the retarded correlator 
of the currents in the medium: 

G,,{x -x\t- t') = -ie{t - t') {[Mx,t), Jj{x',t')]) . (8) 



'Notice that the magnetic permeability fi can lose its significance at relatively small frequencies [21] 
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We use for G the same decomposition in transverse and longitudinal part used in eq. ([3]), 
and using the macroscopic Maxwell equations we obtain: 



where q is the four- dimensional EM coupling. Expanding to second order in k: Gt{u}, k) = 
G'^\u) + k'^G^^\u), we find the electric permittivity and effective magnetic permeability 

as 



where in the second line we have expanded for weak EM coupling q. In summary, the 
response functions, determining the propagation of light in the medium, are given in terms 
of the retarded correlator of transverse current^. 

2 Negative refraction 

As we just discussed, the propagation of electromagnetic waves in a medium is usually 
described in terms of the refractive index n defined as = e/x. In the absence of dissi- 
pation e and fi are real. In this case, there is wave propagation in the material only if 
e/i > 0. The definition of n does not seem sensitive to the simultaneous change of sign of 
both e and /x, since the refractive index is defined as a quadratic equation. However, it 
was understood quite some time ago [H] that changing the sign of e and fj, corresponds to 
changing the branch in the square root, i.e. passing from n = +^JeJl, positive refractive 
index, to n = —^feji, negative refractive index. 

The phase velocity of a wave is defined as Vph = 1/Re(n), and the change of sign 
corresponds to the change of direction of the phase velocity. On the other hand the 
direction of the Poynting vector, and then of the energy flow, is not affected by this 
change of sign. We are in the exotic situation in which the energy flow and the phase 
velocity are opposite. In this case many physical laws change and the material has very 
special properties such as inverse Doppler and Cherenkov effect, inverse Snell's law, the 
medium could work clS cl perfect lens, etc. 

However, the situation with constant and negative e and ^ is non physical [Ij. One 
needs e and to have a dependence on the frequency. In this case of frequency dispersion, 
from general principles we expect e and /i to acquire an imaginary part, so that the medium 
is dissipative. In these media n is itself a complex quantity. Its real part is the index 
of refraction while its imaginary part (usually referred as the extinction coefficient) takes 

^Longitudinal waves can also propagate if there are solutions of the second equation of (|4]), and they 
could also in principle exhibit negative refraction, but in this paper we will only consider the transverse 
modes. 




(9) 




(10) 
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into account the dissipation, and it is always positive in lossy materials. Differently from 
the case of real e and /i, here the refraction can be negative also if Re(e) and Re{^) are 
not simultaneously negative, because of the presence of the imaginary parts. 

For the dissipative case, many different equivalent conditions for negative refraction 
have been worked out [23,i2l]. In [21] (as we review in appendix E]) it was shown that the 
energy flow and the phase velocity are opposite if and only if the index 

riDLioo) = \e{uj)\Re{i2{uj)) + \ 12(00) \Re{e{uj)) (11) 

is negative. In the opposite case we have normal wave propagation. The condition udl < 
is equivalent to require that the refractive index is negative. In this paper we will use udl 
to check if the medium has positive or negative refraction. 



It is natural to wonder for which kind of media tidl < 0. In the following we study 
the behavior of in systems that admit an hydrodynamic description for low frequen- 
cies and long wave vector, and argue that there is a simple condition that implies udl < 
in these systems. We assume the existence of a current J such that its transverse part 
Jt has a diffusive behavior: [dt — W'^)Jt = 0, where V is the diffusion coefficient. This 
implies that the retarded correlator of Jt has a diffusive pole, and it has the form 

with B a real number, up to contact terms and terms that are subdominant in the hy- 
drodynamic limit (i.e. they are of higher order in u, k"^). From generic properties of the 
Green functions we know that: V > (causality implies that the poles can only lie in the 
lower complex half-plane) and B > (Im Gt < 0). If Jt couples to the EM field with 
coupling q, the electric permittivity and magnetic permeability are 

e{uj) = 1 —, i2{uj) = 1 + (13) 

w"^ 00 

From these expressions it is easy to show that n£,L < if < Anq'^B: every system that 
couples to the EM field with a transverse conserved current with retarded correlator as 
in f|T2|) shows negative refraction at low enough frequencies. 

The analysis of the hydrodynamic equations for a relativistic system at finite temper- 
ature and chemical potential |25l|26] shows that the transverse current indeed satisfies a 
diffusion equation; it is also possible to relate the constants B and V to the values of the 
transport coefficients and thermodynamical quantities: 

e is the energy density, p the charge density, P the pressure and rj the shear viscosity. 
Even though it is a simple prediction of hydrodynamics |^, this observation is new, to our 



^ We are grateful to Sean Hartnoll for suggesting us to investigate the generic structure of the corre- 
lators dictated by relativistic hydrodynamics 
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knowledge. It would be interesting to perform a detailed analysis of the EM linear re- 
sponse functions for a generic relativistic hydrodynamic system at finite temperature and 
chemical potential, and give a general claim about negative refraction at low frequencies. 

From now on we concentrate on a specific class of systems: strongly coupled plasmas 
at finite temperature and chemical potential, eventually coupled to a dynamical photon, 
that admit a dual geometric description. They are very interesting for two main reasons: 
it is possible to explicit compute the transport coefficients and the thermodynamical 
quantities, and it is possible to study the transverse current- current response function for 
every frequency, going well beyond the hydrodynamic regime. Even if they are strongly 
coupled systems we can perform the calculations by using the holographic approach. We 
find a range of frequencies where the refractive index is negative, for a generic class of 
these materials. In particular the transverse GT{uj,k) behaves as (fT2!) and the constants 
are given in term of the coefficients in fll4l) . for low frequencies. 

In Figure [T] we show the behavior of the index udl for fixed temperature and for 
different values of the ratio S/T, where S is the chemical potential, and T the temperature 
(fixed to T = 1 for simplicity). 

We observe that, as expected, at low frequencies there is negative refraction, whereas 
for high frequencies there is positive refraction. In the limit of zero chemical potential 
the index tidl is always positive and there is no negative refraction. As the charge is 
switched on, a region of negative refraction appears. This region becomes larger as the 
charge and/or the temperature increase. 

Before giving some details of the system and the computations, it is important to discuss 
the caveat of the dynamical photon that we already outlined in the introduction. It is 
well-known that generically in the holographic correspondence a local U{1) symmetry in 
the bulk corresponds to a global U{1) symmetry in the boundary. As a consequence there 
is no dynamical photon in the dual field theory described by the five dimensional geome- 
try. As a natural way out we could consider our system as a strongly interacting medium 
weakly coupled to the gauge field associated to the global current we are considering. We 
use the holographic correspondence to describe only the medium and we add by hand, in 
a second time, the corrections induced by the presence of a dynamical photon. We would 
like to argue that this is a sensible procedure. The propagator for a dynamical photon in 
a given medium is pTj : 

4:71 

u^e'^iu, k) - P ^^^^ 

where q is the coupling between the current and the photon: qJ^A^, and e'^{oj,k) = 
1 — ^q'^G'ip{uj,k). G\,{uj,k) is the full retarded current-current correlator. It takes into 
account the interactions among the constituents of the medium and the interactions be- 
tween the medium and the dynamical photon. If the medium and the photon are weakly 
coupled the current-current correlator has a natural expansion in power of the coupling 
constant q: C^lu, k) = Gt{uj, k) + q'^G'^\u, k) + ... In this paper we are able to compute 
only the first term of this series: Gt{u}, k). It is the leading term in perturbative expan- 
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Figure 1: for T = 1 and increasing values of the chemical potential. We plot uol as 
a function of the parameter a, implicitly defined by the relation p = The reason 

will become clear in the following sections. 



sion of (|T5l) and we could be confident that we are considering the leading effects for the 
propagation of a dynamical photon weakly coupled to the strongly coupled medium. The 
same philosophy has been applied to the computation of the photo-production rate from 
the quark-gluon plasma in [28] . 

3 The holographic model 

In this section we introduce our model and we explain in which sense our optical systems 
are realized in geometry. 

It is in general quite hard to derive e{u), ^{u) for a certain given material. In this 
paper we take a somehow different road and we compute e^u), niu) for a particular class of 
strongly coupled plasma coupled to an electromagnetic field. The media are characterized 
by two main parameters: the temperature and the chemical potential. 

We use the holographic principle to study these systems: the strongly coupled phase of 
a quantum field theory living in four dimensional Minkowski space time can be described 
in term of classical gravity in a five dimensional Anti de Sitter space (AdSs). The tem- 
perature, and the chemical potential of the field theory are encoded in the temperature 
and charge of a black hole solution in AdSs. 

In this minimal setup our media are described by an energy momentum tensor T^j_i, and 
a conserved U{1) current J^. In the holographic prescription they correspond to a five 
dimensional metric Qmn and a five dimensional vector Am respectively. Their dynamics is 
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described by the Maxwell-Einstein action: 



where 6//^ is the cosmological constant. The ground state of our systems is described by 
the particular solution of the equation of motions of (fTB]l : 

where S is the chemical potential and 

f{u) = {l-u){l + u-au^) (18) 

The solution f ll7p is a charged black hole in AdSs: the event horizon is at u = 1, while 
the AdS boundary is at m = 0. The chemical potential and the temperature of the system 
are encoded into the parameter a and b through the relations 

- « 1 [3~ 

^ = ^Ah;a (19) 



Anb 2 6 V 2 



The propagation of the five dimensional photon defines the retarded current-current 
green function for our (3+l)d media. It is necessary to compute G^\ui) and G^\uj) to 
obtain e{u) and fi{u) defined in the previous section. 

The procedure to compute these functions using holography is well-established [29]. 
To compute the (3+l)d transverse retarded green function Gt{w, k) we need to linearize 
around the solution f|T7j) the equations of motions of f|T6l) for the components of the 5d 
photon Am transverse to the wave vector k. The equations for the photon are coupled 
to the equations for the metric and they decoupled only in the limit of zero charge or 
zero wave vector (a — )■ or — )■ 0). We then need to solve these equations imposing the 
infalling boundary condition at the horizon u = 1. This boundary condition correspond 
to the physical requirement that nothing can escape from the 5d black hole at classical 
level, and it is the retarded prescription for correlators of the (3+l)d medium. From the 
solution for the five dimensional photon we obtain the retarded current- current correlator 
for the media and hence the linear response function e{u) and fi{uj). 

This system of equation was solved numerically in [30] and analytically in the hydro- 
dynamic limit in [3TH33] |^. The hydrodynamic approximation is a series expansion in u 
and k for the equations and solutions. It is valid only for small values of u and k, but it 
offers analytic expressions. In Appendix [B] we give the explicit form of the retarded cor- 
relator at the lowest order in u and fc^. Its leading behavior, necessary for the calculation 



'See also [34] for a previous calculation of these Green functions on 2+1 dimensions. 
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of e{co) and fi{oj), is given by 



Using the equations ( fTOl) 




) 



w8(l + a)2e26 



i 3a/ 



(20) 



e(w) = 1 + 



/tt fi (2-a)2 



1 3a 



) 



6(1 + a)e2 Vw2(l + a) 



/i(w) ~ 1 + g' 




(21) 



The form of Gxiujyk), e(a;), /i(a;), is exactly as in (fT2|) . ( [T3|) . (IT^ . except for an higher 
order term in u that does not influence the conclusion of the previous section. Indeed it 
is interesting to observe that as soon as we turn on the electric charge of the black hole, 
i.e. the medium has a non vanishing chemical potential, the low frequency behavior of 
the real part of e{uj) changes discontinuously to negative value, while the imaginary part 
of /i(a;) acquires a positive contribution. Even if Re(/i(a;)) never becomes negative, the 
dissipation and the fact that Re(e(a;))<0 conspire to give negative refractive index at low 
frequencies 0- In the next section we will go beyond the hydrodynamic regime and we 
will compute correlators and response functions for every frequency u. 

4 Numerics and Plots 

To go beyond the hydrodynamic limit we need to numerically integrate the equations 
of motion for the 5d gauge field and the metric. The functions e{(x)), ii{uj) are obtained 
from the correlators G'^^\u) and G^^\u). We obtain these Green functions expanding 
the equations in the wave vector k up to the order fc^ and numerically integrate the 
resulting equations for every value of uj. As usual the solutions in AdS have IR divergences, 
that correspond to UV divergences in the dual field theory, that we need to subtract. 
The holographic renormalization procedure [32] leaves the freedom to add a term to 
the Green functions of the form c{uj'^ — k"^). The constant c represents a renormalization 
ambiguity that has to be fixed. Our prescription is to require that the electric permittivity 
approaches 1 at large frequencies. This is a natural physical requirement: for infinitely 
rapid variation of the fields the medium does not have the time to adapt and it behaves 
like the vacuum. 

^''We consider a perturbative series expansion in the (3+l)d EM coupling q. The vahdity of the 
perturbative expansion would imply a lower bound on the frequencies u, once all the other physical 
parameters of the systems are fixed. 

The negative quadratic pole in t{uj) is related at least in part to the translation invariance of the 
systems we are considering |20j . Indeed in translationally invariant systems the conductivity (j{uj) has 
an imaginary pole i/oJ, and e(w) = 1 + A'Kia[uj) / lo. However we have checked that the refractive index is 
still negative even if we add a small amount of impurities that break the translation invariance. 
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Im(epsT) 




Figure 2: Imer(co', k) for real u and k, with a = 1 and g = 0.05 



With these Green functions we can study the optical properties of the medium that we 
are describing using gauge/gravity duality and we can check if it shows negative refraction 
of light, once coupled to a dynamical photon. 

Here we plot all the results by fixing T = 1 and by varying the parameters a and 
consequently b. We vary a from a = .5 to a = 1.8, so the chemical potential varies as in 
f fT9|) . from S ^ 3.5 to S ~ 52, i.e. S > T. Nevertheless at higher temperature similar 
results hold. The regime of negative refraction enlarges at higher temperature and reduces 
at lower ones. 

We first have to ensure that the system is in thermodynamical equilibrium. As we 
explain in the Appendix El a system in thermodynamical equilibrium must satisfy the 
requirement Im(er) > 0, and we expect our system to satisfy this requirements since 
is given by a Green function. We have checked this property by plotting this function in 
Figure [2] for real u and k. For this plot we have numerically computed Gt{u}, k) at all 
orders also in k. 

Then we plot the real and imaginary parts of the electric permittivity and of the 
magnetic permeability, in Figure [3] and H] respectively. Notice that Re(/i(a;)) and Re(e(a;)) 
are not simultaneously negative. As we stressed in section [2], also if Re(e)- Re(/i) < the 
energy and phase velocity can be in opposite directions, because of the role played by the 
imaginary parts, and the material has negative light refraction. 

Indeed from Figure [T] we see that the index udl flTT]) is negative at small frequencies 
for different values of the parameter a. In the limit a — )■ we reduce to the uncharged 
case, and the index is positive also at small frequencies. For small charge a, the index 
becomes immediately negative at small frequencies, and then it become positive for larger 
u. By increasing the value of a the region of negative refraction grows. 

A last useful quantity to plot is the ratio Re(n)/Im(?T,), in Figure [5l Indeed this 
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Figure 5: Ratio -Re[n(c(;)]/Im[n(a;)] for different values of a and q = 0.05 



quantity takes into account the ratio between propagation and dissipation . Usually in 
isotropic metamaterials it is common to have strong dissipation in the regime of negative 



refraction 



This implies that 



Re{n) 
Im{n) 



<^ 1. Here we see that this is the case too. 

The situation improves as a; — ?■ 0, where the effects of dissipation get lower. At larger 
frequencies instead the dissipation dominates over the propagation. 

As an aside remark, we have found that at higher frequencies Im(/i(ci;)) becomes neg- 
ative. We stress that it is not a signal of any instability, since eT{uj,k) is always a well 
defined response function for every value of cu, in particular it always has positive imagi- 
nary part. We could interpret this phenomenon as the breakdown of the e-/i approach; it 
would be desirable to have a better understanding of this point. 



Conclusions and future directions 

In the last decade there was a revolution in the practical ability to control the electro- 
magnetic fields in various media and the possibility to study exotic optical phenomena. 
In this scenario negative refraction of light has been shown to be possible in artificial 
metamaterials. One can wonder if there are other physical systems with this property. In 
this paper we concentrated on relativistic field theories at finite temperature and chemical 
potential. 

We guessed that negative refraction is generic in the hydrodynamic limit of these the- 
ories. This behavior is based on the form of the retarded propagator of the transverse cur- 
rent. We confirmed this observation with an explicit example. Indeed we studied the op- 
tical properties of some strongly coupled media from the perspective of the gauge/gravity 
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correspondence. In this case there is a procedure to exactly compute the transport coeffi- 
cients and the response functions in terms of the temperature and the chemical potential, 
and we can obtain quantitative informations about the refractive index. 

We argued that a charged black hole in AdSs is a good laboratory for the study of this 
unusual optical property of three dimensional media. In this background we have found 
the range of frequency at which the index of refraction is negative. Even if we studied 
a specific background, our model is rather generic. Indeed the holographic description 
presented in section [3] is valid for every strongly coupled media with at least a U{1) global 
symmetry. It is important to stress that we are not studying the optical properties of 
the complicate five dimensional curved background. This background is only necessary to 
compute the correlators of the dual strongly coupled 3+1 dimensional field theory, whose 
geometry is fiat. Otherwise non-fiat geometries lead to other non-conventional optics 
phenomena (see for example |37j). 

There are many possible future directions of research. 

On a general ground it would be interesting to find if other amazing phenomena 
coming from the recent optics developments can be usefully realized in string theory, or 
if the ability of string theory to describe some peculiar state of matter could help in optic 
devices. 

The U{1) gauge symmetry in the five dimensional gravity background is typically un- 
derstood as a global symmetry in the dual 3+1 dimensional field theory. We hope that the 
possibility of negative refraction could stimulate the efforts towards a better understand- 
ing of the gauging procedure in AdS/CFT [38]. Positive results in this direction could 
give the possibility to actually compute the ffist quantum electro-dynamical corrections 
to the refractive index. 

As a matter of fact we noted that the dissipation is always stronger than the propa- 
gation at frequencies corresponding to negative values of the index of refraction. It is a 
common property of many isotropic media [36], and some ad-hoc engineering procedure 
have been used to reduce losses in metamaterials. It should be interesting to find an 
analogy of this procedure in the system that we studied. Maybe one can try to lower the 
dissipative effects by studying deformations of the holographic background. 

The media we studied in this paper are homogeneous and isotropic. To make better 
contact with real materials it would be interesting to break homogeneity and isotropy 
and to compute current- current correlators in these less symmetric backgrounds [39tB0]. 
Another possibility is the breakdown of the relativistic symmetry [T6|HT]. 

Negative refraction has been observed also in two dimensional materials [121113], and 
it seems natural to extend our analysis of the optical properties from three dimensional 
media to two dimensional ones at weak coupling and at strong coupling (by using the 
AdS4/CFT3 correspondence). 

Another natural extension is the study of optical properties of material in the su- 
perconducting phase. Indeed negative refraction has been discussed in superconductors 
in [H] , and maybe the same property can be studied in holographic superconductors [13] 
(see also [20] and reference therein). 

We conclude by noticing the analogy between negative refraction in photonics crystals 
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and the motion of the electrons in a bipolar junction in graphene |15]. It can be interesting 
to better understand this analogy and we hope that our work can shed new light on the 
holographic dual of this effect. 
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A Wave propagation and dissipation 

In this appendix we derive the conditions for negative refraction and the consistency 
conditions for the e{u, k) imposed by thermodynamics. 

All along this paper we use the single function e{u, k) to take into account all the inter- 
actions between the medium and the electromagnetic field, the charges and the currents. 
In this setup the Maxwell equations take the following form: 

V-L) = , V-5 = , V AE = -dtB , V AB = dtD (22) 

completed by the constitutive relation, written in the dual Fourier space, D = eE. The 
propagation/dissipation of the electromagnetic field in a dispersive/dissipative medium is 
a bit tricky. 

To determine the propagation of energy in the medium, i.e. the Poynting vector, we 
use the approach to consider the field coupled to an external source that compensates 
the energy losses due to the dissipation. In such a way we can think of the medium as 
almost transparent, and the e is a real function of two real coordinates: u, k. Analyzing 
the propagation of a wave packet we obtain the well-known expression for the Poynting 
vector in the transparency regime for a dispersive medium [21]. To take the dissipation 
of energy into account, we just need to analytically continue the expression to complex 
values for k and e. The final expression for the Poynting vector is: 




(23) 
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The Poynting vector is the energy flow. For transverse waves (123!) reduces to: 

If we expand for small values of /c as in ([5]) we obtain the usual expression for the Poynting 
vector in a dissipative medium in the e — /i formalism: 

S = Re(^^yET\\ (25) 

where we have used the relation k/u = n. There is negative refraction when the phase 
velocity is opposed to the flux of energy: 

Re {n) < 0, Re > (26) 

Working with these relations, one can show that (!26|) is equivalent to n£,L{uj) < 0. This 
is the condition we use in the main text to check if the refractive index is negative, which 
was derived in [23]. 

Let us now consider the constraints imposed by the thermodynamics on the e{u,k). 
The divergence of S gives the rate of change of the local electromagnetic energy in a 
medium. To determine the dissipation we consider a monochromatic wave and we take the 
divergence of fl23l) . The first term of fl23|) gives lie {E*dtD + B*B), while the divergence 
of the second term is zero for a plane wave. The time average of minus the divergence of 
the Poynting vector is the energy inflow from the external source needed to maintain the 
plane wave in the medium. Because in a dissipative medium all the energy of the wave 
will be eventually transformed in heat flow, this quantity is equivalent to the outgoing 
heat flow; for transverse waves one has: 

Q = Im(eT)|^Tr (27) 

Our system, in absence of the external electromagnetic wave perturbation, is in thermal 
equilibrium. Hence the energy flow must be positive: Q > 0. We obtain the consistency 
condition: lm{eT{oo,k)) > 0. Our epsilon functions are derived from retarded Green 
functions in a causal, unitary field theory for which Im {G{u, f0) < 0. For this reason the 
consistency condition for the epsilon is automatically satisfiecO, as shown in Figure [2l 



B Finding the correlators 

In this appendix we describe the calculation of the retarded Green functions that we used 
in the paper [30H33] . We need to solve the equations of motion for the fluctuations of 

^^The inequality holds for k real and uj lying in the upper complex half-plane; again this condition is 
off-shell from the point of view of a photon propagating in the medium but we imagine to work with an 
external field that maintains the electromagnetic wave. 
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the metric and the Maxwell field around the solution (fT7|) . However these equations are 
coupled and it is better to define the new variables 



^^iu) = -hfiu) - -iaBJu) + —B^{u) (with C± = (1 + a) ± J{T+af+^dFi^) 

u u 

(28) 

h^{u) is a metric fiuctuation and Bx{u) = ^^^"-^ , with A^^u) a gauge field fiuctuation. The 
equation of motions for the master variables $± are 

^'±i^) + W + -4r-, - ^V(«)) ^±iu) - ^Mu) = (29) 

u^f{u) uf\u) f[u) 

It is necessary to solve these equations and find the expression for B'^{u) to obtain the 
lorentzian correlators. 

The analytical solution is found at the lowest order in u and k'^. First we must impose 
the infalling boundary condition at the horizon u = 1. This constrains the solution to be 

of the form $±(m) = (l — u) f^-") The functions 0± are found by solving the equations 
at the lowest order. After imposing the boundary conditions for all the fields at m = 0, 
the derivative of Bx{u) at the boundary is given by 

= 2zc.(l + a)-&P + ' Ail + ay + ^ ^ ^^^'^^^^^ ^''^ 

The retarded Green function is obtained from the boundary action 

Srt/ f d^k 

by using the prescription of [29] for the minkowskian correlator, and subtracting the UV 
divergencies: 

^ _ ( 3a {2-aY \ 

e2 V262(22a;(l + a)-6P) 8(1 + a)26y' ^ ' 

To obtain a result valid for every value of uj and fc, we need to implement the above 
procedure at numerical level. 
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